Casimir-Polder potential and transition rate in resonating cylindrical cavities 
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We consider the Casimir-Polder potential of particles placed inside a metallic cylindrical cavity 
at finite temperatures, taking account of thermal non-equilibrium effects. In particular, we study 
how the resonant (thermal non-equilibrium) potential and transition rates can be enhanced by fine- 
tuning the radius of the cavity to match the transition wavelength of the dominant transition of 
the particle. Numerical calculations show that the CP potential acting atoms prepared in low-lying 
Rydberg states can be enhanced beyond 30kHz, which is within the range of observability of modern 
experiments. Because the magnitude of the resonance peaks depend sensitively on the low frequency 
dissipation of the cavity metal, experiments in this set-up could be a critical test of the disputed 
thermal correction to the Casimir force between metal plates. 

PACS numbers: 34.35.-|-a, 12.20.-m, 42.50.Ct, 42.50.Nn 



I. INTRODUCTION 

Casimir-Polder (CP) forces [1] belong to the group of 
of dispersion forces which arise due to the fluctuations of 
the quantized electromagnetic field. They occur between 
polarizable atoms or molecules and metallic or dielectric 
macroscopic bodies and can be intuitively thought of as 
the dipole-dipole force caused by spontaneous and mutu- 
ally correlated polarization of the atom or molecule and 
the matter comprising the body. Under the assumption 
of thermal equilibrium, CP forces have been commonly 
investigated in the linear-response formalism [2-4]. 

Recent theoretical predictions [5] as well as experi- 
mental realizations [6] for CP forces in thermal non- 
equilibrium situations have pointed towards interesting 
effects which arise when an atom at equilibrium with its 
local environment interacts with a body held at a differ- 
ent temperature. In particular, depending on the tem- 
peratures of the macroscopic body and the environment, 
the force can change its character from being attractive 
to repulsive and vice versa. 

Non-equilibrium situations between atom and local en- 
vironment can be investigated by means of normal-mode 
quantum electrodynamics (QED) [7, 8] or macroscopic 
QED in absorbing and dispersing media [9, 10]. In this 
case, thermal excitation and de-excitation processes lead 
to resonant contributions to the force [11] (cf. similar 
findings reported in Ref. [12]). At particle-body separa- 
tions that are larger than the wavelengths associated with 
the dominating atomic transitions (retarded regime) , the 
interaction potential becomes spatially oscillating [13]. 
Similar behaviour has been observed for the transition 
rate of molecules in the past [14, 15]. 

The spatially oscillating non-equilibrium forces on 
ground-state atoms or molecules are proportional to the 
thermal photon number. In order to observe them, it 
is necessary to make use of atomic systems whose inter- 



FIG. 1: Cross-section of the geometry considered: a particle 
in a vacuum-filled circularly cylindrical cavity. 



nal eigenstates exhibit energy separations of order k^T 
or less (T: temperature, fee: Boltzmann constant). We 
have investigated polar molecules with their low-energy 
rotational and vibrational transitions as possible candi- 
dates [13]. However, the large photon numbers obtained 
using molecules with small excitation energies come at 
the cost: Due to the large wavelengths associated with 
such molecular transitions, the retarded regime where 
oscillations might be observed sets in at very large dis- 
tances, typically of the the order of tens to hundreds of 
microns. As the CP potential decays away rapidly from 
the body surface, it is very small at such distances. 

As already discussed in conjunction with excited atoms 
in a zero temperature environment, resonant forces can 
be enhanced in a planar cavity whose width is fine- 
tuned to match the wavelength of the transition [16-22]. 
Note that resonating cavities have been employed exper- 
imentally for enhancement and inhibition of spontaneous 
emission rates for excited systems for a long time (cf. e.g. 
[29-32]). We have previously investigated the potential 
of such a setup to enhance the predicted spatial oscilla- 
tions of the thermal force on ground-state molecules [23] . 



Unfortunately, the cavity enhancement factor turns out 
to scale logarithmically with the cavity Q factor, which 
strongly limits the possibilities of the scheme. As we have 
shown, it is unlikely to achieve more than an order of 
magnitude's enhancement of the force amplitudes, which 
is still insufficient for detection using polar molecules. 

In the present article, we present a scheme by which 
oscillations of the resonant thermal CP force can be 
brought into the measurable regime. This is achieved 
by replacing the planar cavity with a cylindrical cav- 
ity [24] and employing Rydberg atoms rather than polar 
molecules [25]. The geometry is shown in Fig. 1. The 
limited cavity enhancement in a planar setup is due to 
the insufficient, purely one-dimensional confinement of 
the electromagnetic modes. A cylindrical cavity is there- 
fore an obvious candidate for improvement: It confines 
the modes in two out of three spatial dimensions, so that 
a stronger resonant enhancement may be expected. At 
the same time, it is a practical geometry for experimental 
and guiding purposes, allowing particles to travel freely 
along the axial direction (we will only consider particles 
at rest herein, leaving moving atoms [26] in this geometry 
for future study). As a caveat, one has to bear in mind 
that a cylindrical cavity will simultaneously enhance both 
the CP potential and the relaxation rates; thus reducing 
the time scales on which non-equilibrium effects can be 
observed. A similar geometry, i.e., an (anisotropic) par- 
ticle above a spherical hole in a thin metal plate, has 
recently been suggested as a candidate for observing re- 
pulsive Casimir forces [27]. 

Rydberg atoms are ideal candidates for observing reso- 
nant thermal CP potentials. Their valence electrons are 
excited to relatively stable states with very large prin- 
cipal quantum numbers n, typically in the range 30-60. 
The spatial extent of such highly excited atoms is then 
enormous on an atomic scale, exceeding a micron in di- 
ameter, and the transition dipole moments consequently 
orders of magnitude larger than those of ground-state 
atoms or polar molecules. Resonant CP interactions be- 
ing proportional to the respective transition dipole mo- 
ment squared, a strong enhancement follows. Moreover, 
an atom excited to a particular eigenstate with large n is 
necessarily out of thermal equilibrium with its environ- 
ment, and the energy difference to neighbouring states 
is typically small compared to k^T at room tempera- 
ture, hence fulfilling the condition for observing oscillat- 
ing potential. Remarkably, we recently found that the 
same property ensures that the CP potential in a Ryd- 
berg atom is virtually independent of temperature from 
room temperature down to absolute zero [28] . 

Outline of the paper. The structure of the paper is 
as follows. In section II we present the general formal- 
ism of the thermal CP interaction and transition rates in 
a cylindrical cavity, including the necessary Green ten- 
sor. Thereafter, in section III we derive the cylinder radii 
which resonate with the atomic transition frequency. We 
begin with the simplest case of perfectly conducting cav- 
ity and then discuss how the optimal radii deviate from 



the perfect conductor results when realistic metallic per- 
mittivity functions are employed. We provide simple for- 
mulas for the optimal radii when the permittivity of the 
cylinder medium is large. In section IV, we discuss how 
the potential and heating rate enhancements scale with 
the relevant atomic and cavity parameters, with empha- 
sis on the dependence on the cavity permittivity. Fi- 
nally we undertake numerical calculations of resonant 
CP potential and heating rates for two example cases: 
the 32si/2 — >■ 31p3/2 transition of Rydberg Rb in a cylin- 
drical Au cavity at temperature T — 300K, and a ground 
state LiH molecule in a similar cavity. 



II. GENERAL FORMALISM 

We consider an atomic system (below we study the 
cases of Rydberg atoms and polar molecules) with inter- 
nal energy eigenstates ]n), eigenenergies fuvn, transition 
frequencies cvmn = '-^m — ^n and dipole matrix elements 
dmri = (?Ti]d]n), which is prepared in an incoherent su- 
perposition of its energy eigenstates with probabilities 

Pn- 



A. Thermal Casimir-Polder potential 

The thermal CP force on such a system in an envi- 
ronment of uniform temperature T was derived in detail 
in Ref. [11]. As shown, it is conservative in the pertur- 
bative limit, F(r) = — VC/(r), where the associated CP 
potential is given by 



;7(r)=^p„[/„(r). 



(2.1) 



The potential components associated with a given eigen- 
state n splits naturally into a non-resonant contribution 
f/f"' and a resonant one UI,'^^, 



t/„(r) = C/-(r) + C/r(r). 



(2.2) 



The non-resonant potential [/^"^ is due to virtual photons 
and is reminiscent of that obtained by a dilute-gas expan- 
sion of Lifshitz' formula [2]. The resonant contribution 
U!^'^^ is due to absorption and emission of thermal pho- 
tons; it is present because the particle in its ground state 
is out of thermal equilibrium with its environment. 
The non-resonant potential reads 



K'{r) 






|Tr[a„(*e,)-G(i)(r,r,*e,)] (2.3) 



where /io is the free-space permeability, ^^ — 2'KJkBT/h is 
the jth Matsubara frequency, and the prime on the Mat- 
subara sum indicates that the j = term is to be taken 
with half weight. The atomic/molecular polarizability is 
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kn^nk 



Mik^-^kn 



given by 

a„(w) — lini — y 

(2.4) 
and G^^'(r, r',w) is the scattering part of the classical 
Green tensor of the geometry the particle is placed in. 
Note that dyadic multiplication is implied for products 
of vectors without multiplication symbol: AB = A (g) 
B. For an isotropic particle, the non- resonant potential 
components may be simplified to 



K'{^) 



kBT^,e 



^'4a„(zO)TrG«(r,r,ze,) (2.5) 



3=0 



with 



a, 



,(w) = lim ^ > -^ 2~~ 



3h ^-^ LJ'^^^ — LJ^ — ieuj 



(2.6) 



The resonant potential reads 

[/r(r) = /xo^c^Ldnfc • ReG(iHr,r, |wfc„|) • dfc„ 

X {Q{ujk7i)n{L0kn) - 0(cj„fe)[n(w„fc) + 1]}, (2.7) 

where /xq is the free-space permeability, and 9(x) denotes 
the Heaviside step function. The photon number follows 
the Bose-Einstein distribution 
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(2.8) 



For an isotropic particle, the resonant potential compo- 
nents reduce to 

UTi^) = ^5]c.L|d„fe|^TrReG«(r,r,|^fe„|) 
•^ k 

X {<d{uJkn)n{uJkn) - Q{uJnk)[n{uJnk) + !]}• (2.9) 



B. Environment-assisted transition rates 

An atomic system intially prepared in a given energy 
eigenstate \n) and placed in an environment of uniform 

I 



temperature T will undergo transitions to different eigen- 
states due to an absorption and emission of thermal pho- 
tons. As shown in Ref. [41], the total rate 



r„(r) = r(")+rW(r) 



(2.10) 



for transitions out of state \n) consists of a free-space 

part r„ and an environment-inducet 
perturbative limit, these are given as 



part r„ and an environment- induced part r„ . In the 
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-{Q{uJnk)[niuJnk) + 1] 



+ Q{uJkn)n{uJkn)} 



(2.11a) 



and 



rli)(r) = A V 4" d„fe • ImG(i)(r, r, \ukn\) ■ dfc„ 

k 

X {e(a;„fe)[n(cj„fe) + 1] + e(a;fe„)n(cjfc„)}. (2.11b) 



The environment-induced rate simplifies for an isotropic 
particle to 



k 
X {Q{uJnk)[niuJnk) + 1] + Q{iVkn)n{uJkn)}- (2.12) 



C. The Green tensor in a cylindrical cavity 



We are going to study atomic systems placed at posi- 
tion r — {p, (fi, z) inside a circularly cylindrical free-space 
cavity of radius i? in a bulk non-magnetic medium with 
permittivity s = £(w) as shown in Fig. 1. The respec- 
tive scattering Green can be found in Ref. [33] (see also 
Appendix A4.2 of [10]): 



oo oo 



G^^^(r,r,a;) =— / dg E ^ ^ [?'MMMe„^((7)Mo„^(-g) ± rNMNo„^(g)Me„^(-g) ± rMNMo„^(g)Ne„,,(-(j) 

+rNNN=„^(g)No„,^(-q)] (2.13) 



where 77 = -^//j^ — (p- and fc = uj j c. The cyhndrical vector wave functions inside the cavity are [10, 33, 34]: 



Ne„^((7) 



T— Jmi-qp) ^^^rrnpQ -r]JLivP) ^^^nifip 
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(2.14a) 
(2.14b) 



The compact wave vector notation used here imphes 
Ac Be = AgBe + AqBo etc., and the upper (lower) sign 
in (2.13) corresponds to upper (lower) index e,o of the 
vector wave functions. The reflection coefficients r can 
be found from a system of linear equations as described 
in Ref. [33]. For the single- interface cylindrical cavity 
in a bulk medium as considered here, the result for the 
diagonal coefficients Tm = '"mm and r^ = tnn may be 
written as 



rM.N 






■rM.N, 



with 



A + Bd 



, a^M,N, 



and 



A = -m^{kRf{qRf{e-l) 



Bm —x^x [eh^x — [hi] + ehih)xix + hjx^ 



Bd =xlx [eh\x^ - [e 



l)hi~]xix+fxl] 



(2.15) 



(2.16) 



(2.17a) 
, (2.17b) 

Bn ^xlx'^[eh\x^ - [ehi] + hih)xix + hjxj], (2.17c) 

(2.17d) 



rjR, xi = rjiR, 771 = ^/eiujyoyl? 



Here, x 

£ — £(w) is the permittivity of the cylinder and we have 

defined the shorthand quantities 

h^h{x), hi^h{xi), j = j(x), j^=j(xi) (2.18) 

where the reduced Bessel functions denoted h(x) and j(x) 
are 



7/ N Hm [x) d rj(l)f . 



dx 
5(a:)=:^ = AlnJ„(x). 



H'r:,\x) 
^ 



(2.19a) 
(2.19b) 



Note that j is a real function for real arguments whereas 
h is complex for real arguments. Explicit knowledge of 
the off-diagonal reflection coefficients is not required; for 
our purposes, it is sufficient to note that tmn = ''nm (cf. 
Refs. [10, 33] for details). 

Substituting the vector wave functions (2.14a) and 
(2.14b) into Eq. (2.13), making use of tmn — ?'nm and 
the fact that odd functions of q do not contribute to the 
integral, the Green tensor is found to take the diagonal 
form 



G(^'(r,r,( 
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(2.20) 



r 



The trace of the Green tensor required for isotropic 
molecules hence reads 



TrGW(r,r,c.) = — / dqj^' 

"" -^0 m=0 



TM + ^rN 



Tq2p2 



JniiVP) + J'^i^P) 



rN^JiiVP)}- (2.21) 



It should be noted that the terms containing ni?' /if p^ 
exhibit third-order poles at q = k where 77 vanishes. 
The physically correct treatment of this pole consists of 
adding a small imaginary part to the free-space wave 



vector k = lo/c + iS and performing the limit 5 — > 0. 
In this way, the pole is circumvented from below as q 
is integrated along the positive real axis, see Fig. 2(a). 
Correspondingly, 77 is integrated along a part of the real 
axis, around the singularity at the origin and along the 
positive imaginary axis. Fig. 2(c). Recall that in planar 
geometries the Green tensor at real frequencies may be 
separated into a propagating part (q < k) which exhibits 
an oscillating behaviour due to interference of incoming 
and reflected photons and an evanescent part {q > k) 
which falls off monotonously away from a surface [13]. 
In the present case, however, such a distinction cannot 
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(2.23) 



(2.24) 



FIG. 2: (a) and (b); Rotation of the g-integration path by an 
angle 9. (c) and (d): Resulting paths for 77. 



be made due to the pole at the separation point q = k; 
the individual propagating and evanescent parts would 
diverge. 

For numerical purposes, the original integration con- 
tour is unfavourable since it involves subtraction of al- 
most equal numbers, associated with a considerable loss 
of accuracy. Instead, we rotate the contour of the inte- 
gral over the variable q to lie along a line at a small neg- 
ative angle 9, below the real axis as shown in Fig. 2(b). 
As shown in Fig. 2(d), the corresponding path for 77 is 
slightly shifted with respect to the orginal one but still 
contained within the first quadrant: When <C 1, the 
point of closest approach to the pole at 77 = is ap- 
proximately 9{l+i) and the integration path thence ap- 
proaches an asymptote at an angle 7r/2 — 6. 

We will show in Sec. IIIB that all the poles of the 
reflection coefficients lie below the real x axis when |e| < 
00, approaching the real axis from below in the perfect 
conductor limit. There are no other possible singularities 
in the integrands in Eqs. (2.20) and (2.21), so the area 
contained between the original and deformed integration 
contours for 77 is free of poles. As the intergands vanish 
along a path at imaginary infinity connecting the two 
paths, they are equivalent by virtue of Cauchy's theorem 
as long a.s 9 < tt/2; we choose 9 = 0.1 rad. An additional 
benifit of the rotated integration path is the fact that the 
integral becomes exponentially convergent. 

The Green tensor becomes particularly simple on the 
cavity axis, p — where only a few of the functions 
Jm(j]p)^ J'mivp) ^-rc different from zero. We have for 
m = 0,l,2,... 



^'/™M-^0,i,0,0,... (2.22a) 

^™(w)^0,i,0,0,... (2.22b) 

J^(r?p) ^1,0, 0,0,... (2.22c) 



Finally, let us briefly discuss the Green tensor at 
purely imaginary frequencies w = i^ as required for 
the nonresonant potentials (2.3) and (2.5). We have 
77 = i\/ i^ I (? + (p' = iC, so the arguments of the cylindri- 
cal Bessel and Hankel functions appearing in the reflec- 
tion coefficients (2.15)-(2.17) become purely imaginary. 
One has 

Jm(*2/) =i™/m(y), (2.25a) 

j;„(72/) ^i^-H'„Xy\ (2.25b) 

^i''(*2/)-f^"^™+'^i^m(y), (2.25c) 

ll^^}\^y)^\^-^KAy). (2.25d) 

hence the reflection coefficients become (a — M, N) 



r4^0^^{-lr^rA^0■ (2-26) 



The reduced reflection coefficients are found from (2.16) 
with the substitutions x ^ iy = iC,R, Xi — > iyi = iCiR 
with Ci = Ry/e{i(,)£,'^/c'^ + q"^, and j —i't/i,h —i' k/i with 
reduced modified Bessel functions analogous to (2.19): 



Iniy) 



Kn{y) 



Using the reflection coefficients (2.26), the Green ten- 
sor (2.20) at purely imaginary frequencies is seen to be 
purely real and given by [k, — ^/c) 



1 /-oo ^ t- 



m=Q 
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^ll{Cp)fNii). (2.27) 



Correspondingly, its trace reads 



and so the reflection coefficients (2.15) become 



1 t^^ Q^ ( 






Cp' 



■ iUCp) + CiCp) 
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KrniCp) 

IrniCp) ' 

(2.28) 



Note that for purely imaginary frequencies, the (/-inte- 
gration in the Green tensor is unproblematic as rj remains 
finite along the real q-axis. 



III. RESONANT RADII 

Combining the general expressions for the CP potential 
given in Sec. II A with the Green tensor of the cylindrical 
cavity as laid out in Sec. II C, we can explicitly calculate 
the full potential of any particular atom or molecule in a 
given state placed in a cavity of given size and material. 
Such examples, which require a numerical analysis, will 
be given in Sec. V. 

Our main intention is to enhance the resonant poten- 
tial by means of the cavity. To that end, it is worth 
recalling the case of planar cavity: The resonant CP po- 
tential of an atom placed between two plane surfaces is 
enhanced for certain interplate separations [23]. For a 
given transition frequency there exists a series of such 
resonant separations corresponding to integer multiples 
of the atomic transition wavelength. In such a case, the 
transition is (near-) resonant with a standing-wave mode 
of the planar cavity. Mathematically speaking, the en- 
hancement results from a closest matching of the transi- 
tion frequency with a pole of the scattering Green tensor. 

Similarly, for the cylindrical cavity we find that for a 
given transition frequency there exists a series of discrete 
radii such that the transition is near-resonant with one of 
the cavity modes. In this section, we explore analytically 
the structure of these resonances. 



Perfect reflector 



tm 






rN 



H^rk\x) 



(3.3) 



The potential (2.2) diverges, and is thus truly reso- 
nant, when at least one of these coefficients has a pole at 
(7 = 0. Since x = rjR = kR for h = 0, it is clear that reso- 
nances occur in the perfectly reflecting limit when either 
Jjn{kR) = or J^{kR) = 0, that is, when the radius 
equals one of the radii given as 



i?. 



■mj 



-J 
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m] 



(3.4) 



where jmj and j' ■ are the j th zero of Jm {x) and J^j {x) , 
respectively (only zeros jmj > are considered) . 

For each mode m and each polarisation, there is hence 
a number of possible radii Rmj leading to a resonant 
enhancement of the CP potential. The strongest reso- 
nance is that corresponding to the smallest resonant ra- 
dius, which is R'li, corresponding to the first zero of J(, 
j[i « 1.8411838. If the dominant transition \n) — > |fc) is 
a downward one, i.e. a;j,„ < 0, this resonance corresponds 
to a potential minimum which can act as a guiding po- 
tential. For an upward transition (such as will be the 
case e.g. for a ground state molecule), the strongest reso- 
nance which corresponds to a potential minimum in the 
cylinder center is i? = Rn — R'qi, which happens to be 
a double resonance since ju = Jq^ sa 3.8317060. We will 
analyse these examples further in section V. 



B. Good conductor 

In reality, any metal has a finite conductivity so that 
|e| < oo. As we will now show, this results in a shifting 
of the values of kR which give poles at q ~ away from 
the zeros of the Bessel functions. As a consequence, the 
optimal radii for enhancing the CP potential or transition 
rates deviate from their perfect-conductor values as given 
by Eq. (3.4). We derive approximate formulae for the 
new optimal radii, valid for good conductors. 

In the following we consider a metal described by the 
Drude model, 



Let us begin with the limit |e| — )■ oo of a perfectly 
conducting cavity. For large |e|, the leading order terms 
of the coefficients A and Bm.n.d in Eqs. (2.17) read 
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m^{kRf{qRfe^ + 0{e) (3.1a) 

{kRfi-qRfWhe^'^ + 0{e^) (3.1b) 
{kRfi-qRfhije^'^ + 0{e^). (3.1c) 



Hence, in the limit \e\ 



TM 



> oo, Eq. (2.16) simplifies to 
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(3.2) 



e{^) 



1 



uj{lj + i"f) 
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w + 7 cj o;^ + 7^ 



(3.5) 



where cjp is the plasma frequency and 7 is the relaxation 
frequency. For good conductors one has cjp ^ 7, and 
we assume this is true in the following. Since we are 
considering molecules, we will restrict our interest to low 
frequencies and therefore also assume a; <C Wp. Under 
these assumptions the following is true 



Ree < 0; Ime > 0; |Rce|,Ime>l. 



(3.6) 



Let us consider the reflection coefficients tm^n as given 
by Eqs. (2.15)-(2.17) once more. We are looking for the 



complex resonant value of kR which corresponds to a 
zero of the denominator of tm.n when q = 0. A van- 
ishes quadratically for small q, so we can set A = in 
the following and consider the coefficients Bm.n.d- The 
sought value of kR is hence the solution of the equation 
-Bd|(2=o = [42], which may be written as 



h{y/ekR) %kR) 



1 



~]{kR) h{^kR) 



(3.7) 



The two roots of this second-order equation in ~]{kR) give 
all the resonances since Br, is the dcnomiator of both 
reduced reflection coefficients for q = 0. The perfect 
conductor limit jej — > oo is easily recovered from this 
equation, in which case the solutions ~]{kR) — oo and 
]{kR) = are just the solutions (3.4), i.e., kR ~ jmn and 

KJX — Jfyin- 

Approximate solutions to (3.7) when |e| is large but 
finite are straightforward to find. We write 



kR = j 



(') 

mj 
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(3.8) 



where 5^'' are small complex numbers. The prime cor- 
responds to solutions close to a zero of J^. Solving 
Eq. (3.7) to leading order in e^^ then gives 






S' = 6l 



mj 



^ 'Jm [Jmj ) 1 



(3.9a) 
(3.9b) 



note that h{z) - i + 0{z-^), \z\ -^ oo (cf. e.g. [40] 
§9.2). The fraction JmWmj)/JmU'mj) is a real and neg- 
ative number of order unity which tends asymptotically 
to —1 for large arguments. For a good conductor (3.6), e 
is in the second quadrant of the complex plane, so ^/s is 
in the first quadrant for an absorbing medium, and i/y/e 
is in the first complex quadrant. The shifts S and 5' to 
the poles at g = lie in the third quadrant. 

Similarly, one can show that all poles of the integrand 
of the scattering Green tensor of a well conducting surface 
are displaced from the Bessel zeros rjR — jmj and rjR = 
j'mj by small quantities which lie in the third complex 
quadrant. This means that the poles of the refiection 
coefficients all lie in the lower half of the complex -q plane 
when |e| < oo, justifying the rotation of g-integral path 
mentioned in Sec. II C and shown in Fig. 2. 

For real frequencies, the cavity radius can never be cho- 
sen such that it lies exactly on one of the complex- valued 
resonances. Instead, we will derive optimal real radii 
close to the resonances that maximise the real or imag- 
inary parts of the Green tensor. As these optimal radii 
turn out to be different for the real vs imaginary parts, 
we have to distinguish between radii which maximise the 
resonant potential (2.7) and those which maximise the 
heating rate (2.11b). 



C. Optimal radii for enhancing the potential 

Consider first a resonance associated with a pole of rj^. 
We represent the pole in the form (3.8), and let q and 5 
be small, so that 



X^-qRKi jnj + S 



■ Jnj 



k"^ ■ 



Keeping only leading orders in the small quantities S, q^ 
and 1/y/e, i.e., using Eqs. (3.1a)-(3.1c) with 



Bat - - {kRfirjRfhije'^^^ + {kRY{r]Rffe 



^2^2 



0{e' 



3/2\ 



(3.10) 



the reflection coefficient as given by (2.15)-(2.17) be- 
comes 



TN 



l^H^rk\x) 






JLUmj) q^Jm,/{2k^)-C 



with 



C^6 + t 



(3.11) 



(3.12) 



We have used that hi k, i because of its argument being 
large with positive imaginary part (see the asymptotic 

expansions, section 9.2 of [40]) and noted that Hm {x) « 
iYn{x). 

Now we note that in the integrand of the integral 
(2.20), the term which resonates is the last one, which 
does not have a prefactor q'^. For this term there is no 
other resonating structures in the integrand, and we can 
simply conclude that close to a strong resonance 
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U^°^{r) oc ReG(r,r,a;) oc Im / dqvN 

Jo 



oc Re 



_^ oc Re^— oc Im^- , (3.13) 



where we have made the substitution q — >■ ^2/jmjq/k at 
the beginning of the second line. Explicitly, we have 



Im 



icr 
-1 



1/2 . 

sm 



1 ImC 
■k arctan 

2 ReC 



ReC 



V2^ReY+In?C 



v/Re^C + WC 



(3.14) 



where ReC = 5 - Imje^^/^} ^nd ImC = Reje^^/^}. Dif- 
ferentiating this result with respect to S, we find that the 
maximum lies at 



ReC= -ImC/V3, 



(3.15) 
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FIG. 3: The potential of Rydberg Rb in its 32s state at the 
centre (p = 0) of a An cavity, as a function of radius close 
to the resonance corresponding to the first zero joi ~ 2.4 of 
Jo. The cavity resonantly enhances the contributino from the 
32si/2 — ^ 31p3/2 transition. 7?oi = 799.9ii.m as given by Eq. 
(3.4) is the corresponding perfect conductor resonant radius. 



is much greater than unity as assumed above, the shift of 
the resonant radius away from the perfect-condutor re- 
sult is not negligible. We see in Fig. 3 that the potential 
at the optimal radius is about a factor 2.5 greater than its 
value for kRoi — jqi. The resonant radius is given with 
excellent approximation by Eq. (3.17) which for this ex- 
ample predicts the maximum at i? — i?oi ~ — 187nm. 

We note furthermore that the width of the radius res- 
onances is in the order of 500nm, which is expected to 
be well within the accuracy obtainable for production of 
pipes with diameters in the order of hundreds of microns. 
It is also much wider than surface roughness amplitudes 
of good metal surfaces, indicating that the associated di- 
minishing of the CP-potential enhancement is not ex- 
pected to be important. The narrowness of the peaks are 
thus in the order of one three thousandth of the cylinder 
radius, and we do not expect observation and utilisation 
of the resonant behaviour to be hampered by issues of 
production accuracy. 



so that 



D. Optimal radii for enhancing transition rates 



5 « Im{£-i/2| 



V3 



Re{£-i/2| 



(3.16) 



Note that d is independent of m and j to leading order 
in the large parameter e. The resonant radii associated 
with resonances of r^ are hence given by 

kR « j,„, + Im{£-i/2} - i=Re{e-i/2}_ (3.17) 

For resonances due to rjv/, we present the poles of 
J!m{x) in the form (3.8). A virtually identical procedure 
then leads to 



6' = 51 



mj 



^m (Jmj ) 
'^ m\Jmj > 



(3.18) 



Consequently, the respective resonant radii read 



kR' 



Jmj 



'^rniJm-j) 



1 



J"(/ . y^-^i^-'^'} ' ^^-i^-'^'} 

•^rn Umj ) \ V J 



(3.19) 

As a numerical example, take Rb in its 32s Rydberg 
state whose dominating \n) — 32si/2 -^ \k) = 31p3/2 
transition has frequency cukn ~ —9.013 x 10^^ rad/s and a 
cylinder made of Au using Wp w 1.4 x 10^^, 7 w 5.4x 10^^^ 
[35]. In this case, the shifts of the maxima away from the 
Bessel zero for for resonances due to r^ are found to be 
S « -0.00056. 

In Fig. 3, we show the potential as a function of ra- 
dius close to the smallest resonance at A:i?oi ~ joi, res- 
onating with the smallest downward transition to 31p. 
More details on the specifics of Casimir-Polder poten- 
tials on Rydberg atoms are found in Ref. [25] and sum- 
marised in Sec. V below. It is interesting to note that 
even though gold is a good conductor whose permittivity 



Optimal radii for resonantly enhancing transition rates 
in a conducting cavity can be derived in close analogy 
to the previous section. We again start with resonances 
of r^v as approximated by Eq. (3.11). The transition 
rates (2.11b) close to a resonance are found just as in 
(3.13) to be proportional to 



r£i)(r) ex ImG(r,r,a;) oc ReW- 
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where 
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(3.21) 



The maximum of the above function is again found by 
differentiation w.r.t. S. It now lies at Re^ = ImC/\/3, so 
that 
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Im{£-V2} + -^Re{e"V2} 



(3.22) 



is again independent of m and j to leading order in e. 
The resonant radii associated with resonances of rjy for 
enhancing transition rates 



kR. 
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Jmj 



lm{e-^/^} + ^Re{e-^/^} 



(3.23) 



are thus different from the corresponding radii (3.17) for 
enhancing the CP potential. 

As for the potential, the resonances due to tm are 
found to be maximal at 
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FIG. 4: Same set-up as Fig. 3, but for the cavity-assisted 
transition rate from state 32s to 32p3/2- 
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(3.25) 

The resonant enhancement of the total transition rate 
out of the 32s state of Rydbcrg Rb placed at the cen- 
tre of a Au cavity is shown in Fig. 4, again for the joi 
resonance of vn- The true optinal radius is seen to be 
well aproximated by Eq. (3.23), it is smaller than the 
perfect-conductor value by 0.004%. The optimal radii 
for enhancing potential vs rate thus differ noteably by 
about 192 nm. 



IV. GENERAL SCALING PROPERTIES 

In this section we disucss the scaling properties of res- 
onant thermal CP forces, i.e., their dependences on the 
relevant molecular, material, thermal and geometric pa- 
rameters. These were discussed in detail for the case 
of a planar cavity in Ref. [23]; and many of the re- 
sults remain valid also in the cylindrical geometry: The 
resonant potential corresponding to a dipole transition 
from state n to state k is proportional to the absolute 
square of the transition dipole moment, |dfc„p. Like- 
wise, the resonant potential scales with temperature just 
as in the planar case; proportional to the photon num- 
ber n(|cLifc„|). For temperatures larger than the transition 
frequency the temperature dependence is thus approxi- 
mately linear, {/""^^{p) ex T, kuT :s> huJkn- In the opposite 
limit the resonant potential is exponentially suppressed: 
U'^^ip) ex eM-f^kn/kBT), fceT < hcokn- 

On the contrary, the scaling of the potential with the 
permittivity of the cavity walls depends strongly on the 
geometry. For the cylindrical cavity, this scaling follows 
immediately from the discussion of the resonant radii in 
Sec. Ill C. As shown by Eq. (3.13), the resonant potential 
is proportional to \/l/C- This function (3.14) has its 
maximum at Re<^ = — Im(^/v3, where it takes the value 
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FIG. 5: The dependence of the peak potential at resonance 
on 6* = arg(e). The solid graph is [/"^"^ calculated for a 32s 
Rydberg atom whose 32s — )■ 31p transition is in resonance 
with the joi zero of Jo (a;), varying the phase of e throughout 
the physical range but keeping |ej constant and equal to its 
value using Eq. (3.5). The dashed graph is the theoretical 
curve based on the approximate scaling of Eq. (4.1). Both 
functions have been normalized by their values at arge — ■k/2. 



V^/(2V2EIC) with ImC = Reje-^/^}^ To leading order 
in the large quantity e, the potential for near-resonant 
radius hence scales as 



Urir) (X 



11/4 



cos[^ arg(£)] 



(4.1) 



Note that this scaling is an improvement with respect to 
the enhancement achievable in a planar cavity for which 
we found the potential to be proportional to Ine [23]. 
The resonant enhancement in the above approximation 
diverges as arg(e) — >■ tt. A full numerical investigation 
shows that the potential in fact remains finite in this 
limit, cf. Fig. 5. In the displayed example for the arg(£)- 
dependence of the potential of Rydberg Rb in its 32s 
state, the value of the potential changes by more than 
a factor of three as the phase of e is varied from 7r/2 
(purely imaginary permittivity) to tt (purely real per- 
mittivity). This implies that decreasing the dissipation 
rate of the cavity material can increase the enhancement 
significantly. Note that with 7 w 5.4 x 10^"^ [35], the ac- 
tual phase of Au at the relevant transition frequency is 
close to 7r/2. 

The shown dependence of the resonant potential on 
the phase of the permittivity is also interesting in light 
of the thermal anomaly of the Casimir effect for metals 
[36, 37]. This dispute centers in an essential way on the 
description of the dissipation of the metal: employing 
the standard Drude model (3.5) with measured optical 
dissipation data gives a different prediction of the force 
at high temperature than that using a non-dissipative 
plasma model in which one sets 7 = at the outset. Ex- 
periments appear to favour the latter approach [37, 38]. 
The cavity enhancement of the CP potential as a related 
quantum vacuum effect is good system for investigating 
the thermal anomaly further: As Eq. (4.1) shows, the po- 
tential at resonance using a plasma model predicts a sig- 
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nal more than three times that calculated for the Drude 
model (note that |e| is also larger using a plasma model, 
adding to the relative difference in prediction). 

Finally, let us consider the scaling of the potential with 
the transition frequency. The frequency infuences the po- 
tential (2.7) in three ways: Firstly, there is a prefactor 
a;^^n(|ajfcTi|), independently of the particular cavity ge- 
ometry considered. Secondly, there is an additional fac- 
tor of ujnk because the magnitude of oscillations fall of 
as Jm^Jm ^ l/(^ ~ P) away from the cavity walls and 
the resonant cavity radius is in turn proportional to w^n. 
This geometric frequency-dependence is closely similar 
to that found for a planar cavity [23] ; it is a consequence 
of the general scaling law for the Green tensor as estab- 
lished in Ref. [39]. Thirdly, the e-scaling introduces an 
additional frequency dependence. For LOkn ^ 7, Eq. (4.1) 
leads to a yY/uikn scaling. Combining these three ef- 
fects, the peak resonant potentials scale as a;^;J^n(|a;fc„|). 
For comparison, the scaling w^^n(|wfe„|) for a planar 
cavity is slightly stronger due to the much weaker e- 
dependence. 

For atoms in highly excited Rydberg states (cf. Sec. V 
below), the CP potential is dominated by transitions to 
neighbouring states. For these transitions, the frequen- 
cies and dipole matrix elements depend in a simple way 
on the principal quantum number n of the initial atomic 
Rydberg state: For sufhcienty large n, the transition fre- 
quencies can be given as Wfe„:=2Ry/(fi,n'^) (Ry, Rydberg 
energy) while the dipole moments scale as n^. Combin- 
ing this with the dependencies discussed above, one finds 
that the CP potential of a Rydberg atom scales as 



-4.25 



n{LOkn 



-1.25 



n » TIT] 
exp [— (riT/n)^] n <^ nx. 



(4.2) 



where we have introduced a characteristic thermal prin- 
cipal quantum number 



TlT 



2RyY 

kBTj 



(4.3) 



The maximum potential is found for states with a princi- 
pal quantum number around n ~ ut- At room temper- 
ature, riT ~ 10.2, hence the maximum potential is found 
for principal quantum numbers below the Rydberg range. 
For this reason we have chosen a low Rydberg state, 32s, 
for our numerical examples above and in the following. 



V. NUMERICAL RESULTS 

We are now going to present numerical studies of the 
CP potential inside a cylindrical cavity based on the ex- 
act formulas presented in Sec. II. Due to the complexity 
of the formulaic apparatus for the case of the cylinder, 
it is necessary to first ascertain the correctness of the 
numerical calculations. As a numerical benchmark, we 



verified that for positions sufficiently close to the cylin- 
der wall {R — p <^ R) the potential tends asymptotically 
to that outside a half-space (cf. e.g. [13]). 

We arc interested in an observable resonant enhance- 
ment of the potential predicted for the radii derived in 
Sec. III. As clear from the previous Sec. IV, this requires 
a cavity made of out a good conductor and an atomic 
system whose transition frequencies lie close to the peak 
of the thermal spectrum lo^.^ n{\ujkn\) at room temper- 
ature and whose respective dipole matrix elements are 
large. Rydberg atoms with their enormous matrix el- 
ements and relatively small transition frequencies fulfill 
both of these requirements. 

In practice, Rydberg states with principal quantum 
number n in the range 30 to 50 can readily be prepared 
using standard Rydberg lasers. As shown at the end of 
Sec. IV, the optimal choice of n for measuring the reso- 
nant enhancement is smaller than the standard Rydberg 
regime, around n = 10: While transition dipole moments 
increase with higher n, transition frequencies of the dom- 
inant transition decrease further away from the optimal 
frequency value ^ k^T, resulting together in a low op- 
timum n. As a compromise, we use the value n = 32 
for our numerical examples, being a level in the lower 
part of the Rydberg spectrum while still being readily 
available with standard equipment. In all calculation the 
temperature was 300K, although in all the cases consid- 
ered herein there is virtually no discernable temperature 
dependence (cf. also [28]). 

The resonant potential is calculated only for the domi- 
nating transition to 31p3/2 - transitions to higher and 
lower levels contribute significantly only in the non- 
retarded regime close to the cylinder walls which we are 
not interested in in the present investigation (but see [25] 
for details). Results are shown in Fig. 6 for the first eight 
resonances of the coefficient rj^ (corresponding to zeros 
of J„i{x)), and tm (corresponding to zeros of J^(a;)). 
Note how some radii resonate both with tm and tat since 
J'^{x) = -Jr{x). 

For the considered case of a resonantly enhanced down- 
ward transititon, we observe that a potential at the res- 
onant radius corresponding to j^ ■ has j local minima. 
Potentials for jmj have minima on the cylinder axis for 
even m and maxima for odd m. For the j^ resonances, 
the situation is reversed. The double resonances are dom- 
inated by the jmj contribution, so the first of the two 
rules applies. For a given m, the maximum potential 
depth decreases with j. 

One might expect that smaller radii give the deeper 
potential minima, based on the fact that the amplitude 
of oscillations decrease away from a boundary (outside 
a planar half-space the oscillation amplitude decreases 
proportional to inverse distance [13]). Unlike the planar 
case [23] , however, this is not true in general. The lowest 
Bessel zeros are j[^ « 1.8412, joi ~ 2.4048, j'^i « 3.0542 
and jgj^ — ill K, 3.8317. The two first on this list indeed 
correspond to the two largest potential extrema in the 
same order, but J21 in fact represents the shallowest of 
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FIG. 6: (Left panel) Casimir-Polder potential at the eight smallest resonant radii acting on isotropic Rb in the state 32s, whose 
dominating transition to 31p3/2 resonantes with various modes of a gold cavity described by Eq. (3.5). The radii are chosen 
according to Eq. (3.17) (joi, jo2, jii = joi, J12 = J02, and J21) and Eq. (3.19) (jn, j'x2, j'li)- (Right panel) The unenhanced 
potential outside a plane half-space. For comparison, the oscillating part of the shift due to the transition to 31p3/2 alone is 
shown as a dotted line. 
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FIG. 7: (Left panel) Enhanced transition rates at the eight smallest resonant radii acting on isotropic Rb in the state 32s, whose 
dominating transition to 31p3/2 resonates with various modes of in a gold cavity described by Eq. (3.5). The radii are chosen 
according to Eq. (3.23) (joi, J02, jii = joi. J12 = J02, and J21) and Eq. (3.25) (jn, j'^2, j'zi)- (Right panel) The unenhanced 
shift outside a plane half-space. The free space transition rate is shown as a dotted line. 



the eight enhanced potentials considered in Fig. 6. 

The corresponding enhancement of the transition rate 
from 32s to 31p3/2 is shown in Fig. 7 where the cylinder 
radii are picked according to Eq. (3.23) and Eq. (3.25) 
as appropriate. The figure thus shows the maximally 
enhanced transition rate between these two levels. The 
transition rate is increased by about a factor 7.4 for the 
smallest cavity corresponding to j[i compared to the 
transition rate in free space, which for the 32s state is 
about 65 kHz. In comparison, the transition rate outside 
a half-space oscillates about its free-space value within 



variations of about 10 kHz in the far zone, hence showing 
that the cavity can enhance the oscillating contribution 
to the rate alone by about a factor 50. The lifetime of the 
initial state against spontaneous and stimulated decay is 
thus reduced from its free-space value of 50[as to 1.4|j,s 
inside the cylindrical cavity for the j'n resonance. Recall 
from Sec. HI that the potentials and transition rates peak 
at different radii. In Fig. 7, we have chosen the optimal 
radii for enhancing the rates. By contrast, the transition 
rates at the optimal radii for enhancing potentials (not 
displayed) are smaller by about a factor 2. 
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FIG. 8: Comparison of enhanced and non-enhanced heating 
rate of ground state LiH. The heating rate outside a half-space 
[41] shown together with the two highest cyhnder resonances 
(radii from (3.23) and (3.25)) along a central cylinder cross- 
section. 



For comparison, we consider the case of the polar 
molecule LiH. Polar molecules also exhibit transitions 
in the frequency regime ideal for enhancement, but with 
the respective dipole matrix elements being much smaller 
than those of Rydberg atoms. LiH has a lowest tran- 
sition frequency (rotational) of wio = 2.79 x lO^^rad/s 
and corresponding transition dipole moment |diop — 
3.85 X lO^^^C^m^. As a result of the smaller dipole 
moment, the CP potential of polar molecules is much 
smaller than that of Rydberg atoms. In spite of a rela- 
tive enhancement of about a factor 50-100 in a cylindri- 
cal cavity, the achievable potential depth is still less than 
IHz, which is insufficient for moelcular guiding purposes. 

While a cavity-enhancement of two orders of magni- 
tude is modest for the CP potential of a polar molecule, 
it is of interest when regarding the ground-state heating 
rates. As Fig. 8 shows, the transition to a particular state 
(in this case the lowest rotational state) can be enhanced 
by a factor of about 30, considerably reducing the life- 
time of the rotational ground-state against heating from 
2.1 s to 0.064 s. 



oscillating behaviour is observed for the transition rates 
between eigenstates. 

Formulas for calculating appropriate radii for maximal 
enhancement of potential and rates have been derived. 
For a perfectly conducting cylinder, the resonant radii 
are exactly given as the zeros of the Bessel function Jm 
and its derivatives divided by the wave number of the res- 
onating transition. The optimal radii are slightly shifted 
when the cavity is not perfectly conducting, and simple 
expressions for the correction have been derived for good 
conductors. The corrected optimal radii for enhancing 
potential vs transition rates are now slightly different. 

We have shown how the cavity enhancement scales 
with the relevant paramaters of the set-up. In particular, 
we have paid attention to the dependence of the enhance- 
ment on the permittivity £(a;fe„) of the cavity material 
We have shown that the potential at resonance scales as 
[/ oc |£(a;fe„)p^, which is a noticeable improvement over 
the planar geometry, for which the scaling was found to 
be logarithmic, U oc ln£(a;A;„) [23]. A strong dependence 
is also found on the complex phase arge(a;fe„) of the per- 
mittivity. This is interesting in light of the controversy 
surrounding the temperature correction to the Casimir 
force between metals, for which the complex phase of e 
for the metal in question is of the essence. A precision 
experiment of the potential enhancement in a cylindrical 
cavity could be a critical experiment in this respect. 

The cases of a Rydberg atom and a ground state LiH 
molecule have been studied numerically, both of which 
are of experimental and technological interest. We have 
found that the deepest potential minima for Rydberg 
atoms can be obtained for quantum numbers in the lower 
end of the Rydberg regime. With the smallest cavity en- 
hancement (corresponding to the first zero of J{(a;)), a 
guiding potential depth in excess of 30 kHz is obtain- 
able, which is within the region of observability of mod- 
ern experiments. The enhancement factor obtained is ap- 
proximately 50, an order of magnitude better than what 
we obtained for a planar gold cavity [23]. For the po- 
lar molecule LiH, the cavity enhancement was found to 
be insufficient to bring the potential into the observable 
regime. Instead, a considerable enhancement of ground- 
state heating rates can be achieved. 



VI. CONCLUSIONS 



We have presented the detailed theory for a particle 
(atom or molecule) in an eigenstate inside a cylindrical 
cavity carved out of a homogeneous and non-magnetic 
material. A particle out of thermal equilibrium with its 
environment is subject to spatially oscillating Casimir- 
Polder forces near surfaces, and we have focused partic- 
ularly on a scheme to enhance the oscillating force com- 
ponents by fine-tuning the cavity radius to resonate with 
the particle's internal transition wavelength. A similar 
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